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We analyze the three-body charmless decay -B ± — > K ± K ± K^ using a sample of 226.0 ± 2.5 
million BB pairs collected by the BABAR detector. We measure the total branching fraction and CP 
asymmetry to be B = (35.2±0.9± 1.6) x 10" 6 and A CP = (-1.7±2.6± 1.5)%. We fit the Dalitz plot 
distribution using an isobar model and measure the magnitudes and phases of the decay coefficients. 
We find no evidence of CP violation for the individual components of the isobar model. The decay 
dynamics is dominated by the K + K~ S-wave, for which we perform a partial-wave analysis in the 
region m(K + K~) < 2GeV/c 2 . Significant production of the /o(980) resonance, and of a spin zero 
state near 1.55GeV/c 2 are required in the isobar model description of the data. The partial-wave 
analysis supports this observation. 

PACS numbers: 13.25.Hw, 11.30.Er, 11.80.Et 



I. INTRODUCTION 

Charmless decays of B mesons provide a rich labora- 
tory for studying different aspects of weak and strong 
interactions. With recent theoretical progress in under- 
standing the strong interaction effects, specific predic- 
tions for two-body pseudoscalar-pseudoscalar, B — > PP, 
and pseudoscalar-vector, B — > PV, branching fractions 
and asymmetries are available Q, Q, y, Q and global fits 
to experimental data have been performed @ • Im- 
proved experimental measurements of a comprehensive 
set of charmless B decays coupled with further theoreti- 
cal progress hold the potential to provide significant con- 
straints on the CKM matrix parameters and to discover 
hints of physics beyond the Standard Model in penguin- 
mediated b — > s transitions. 

We analyze the decay B ± — * K ± K ± K T , dominated 
by the b — > s penguin-loop transition, using 226.0 ± 2.5 
million BB pairs collected by the BABAR detectorH at 
the SLAC PEP-II asymmetric-energy B factory Q op- 
crating at the T(45) resonance. BABAR has previously 
measured the total branching fraction and asymmetry 
in this mode [j| and the two-body branching fractions 
B ± -> #±0(1020), B ± -> A compre- 

hensive Dalitz plot analysis of B ± — > K ± K ± K T has 
been published by the Belle collaboration |l2T |. 



II. EVENT SELECTION 

We consider events with at least four reliably recon- 
structed charged-particle tracks consistent with having 
originated from the interaction point. All three tracks 
forming a B ± — > A ± A ± A =F decay candidate are re- 
quired to be consistent with a kaon hypothesis using a 
particle identification algorithm that has an average effi- 
ciency of 94% within the acceptance of the detector and 
an average pion-to-kaon misidentification probability of 
6%. 



*Also at Laboratoire de Physique Corpusculaire, Clermont- 
Ferrand, France 

tAlso with Universita di Perugia, Dipartimento di Fisica, Perugia, 
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J Also with Universita della Basilicata, Potenza, Italy 



We use two kinematic variables to identify the sig- 
nal. The first is AE = E — y / so/2, the difference be- 
tween the reconstructed B candidate energy and half 
the energy of the e + e~ initial state, both in the e + e~ 
center-of-mass (CM) frame. For signal events the AE 
distribution peaks near zero with a resolution of 21 MeV. 
We require the candidates to have \AE\ < 40 MeV. 
The second variable is the energy-substituted mass 
toes = V( s o/ 2 + Po -Pb) 2 /EI -p|, where p B is the 
momentum of the B candidate and (Eq,pq) is the four- 
momentum of the e + e~ initial state, both in the labora- 
tory frame. For signal events the toes distribution peaks 
near the B mass with a resolution of 2.6 MeV/ c 2 . We de- 
fine a signal region (SR) with to E s € (5.27, 5.29) GeV/c 2 
and a sideband (SB) with to es e (5.20, 5.25) GeV/c 2 . 

The dominant background is due to events from light- 
quark or charm continuum production, e + e~ — > qq, 
whose jet-like event topology is different from the more 
spherical B decays. We suppress this continuum back- 
ground by requiring the absolute value of the cosine of 
the angle between the thrust axes of the B candidate and 
the rest of the event in the CM frame to be smaller than 
0.95. Further suppression is achieved using a neural net- 
work with four inputs computed in the CM frame: the 
cosine of the angle between the direction of the B can- 
didate and the beam direction; the absolute value of the 
cosine of the angle between the candidate thrust axis and 
the beam direction; and momentum-weighted sums over 
tracks and neutral clusters not belonging to the candi- 
date, J2iPi an d J2iPi cos2 ®h where pt is the track mo- 
mentum and 0i is the angle between the track momentum 
direction and the candidate thrust axis. 

FigureHsb-Ows the toes distribution of the 9870 events 
thus selected. The histogram is fitted with a sum of a 
Gaussian distribution and a background function having 
a probability density, P(x) cx xV 1 — x 2 exp (— £(1 — x 2 )), 
where x = 2toes / v^o an d £ is a shape parameter |l3| . 
The binned maximum likelihood fit gives \ 2 ~ 104 f° r 
100 bins and £ = 21.1 ± 1.6. The ratio of the integrals 
of the background function over the signal region and 
the sideband yields an extrapolation coefficient R q q = 
0.231 ± 0.007. The expected number of qq background 
events in the signal region is — R q q(n SB — n ^g) — 
972 ± 34, where n SB = 4659 is the number of events in 
the sideband from which we subtract the number of non- 
signal BB background events nzP= = 431 ± 19, estimated 
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FIG. 1: The Toes distribution of the 9870 selected events, 
shown as the data points with statistical errors. The solid 
histogram shows a fit with a sum of a Gaussian distribution 
(too = 5.2797 ± 0.0007 GeV/c 2 , a = 2.64 ± 0.07 MeV/c 2 , N = 
2394 ±63) and a background function [l3j. shown as a dashed 
histogram. The shaded regions correspond to the signal region 
and the toes sideband defined in the text. 



using a large number of simulated exclusive B decays. 
The expected number of BB background events in the 
signal region is n s ^- = 276±20, with B^ — > DK ± decays 
giving the largest contribution. 

We use a kinematic fit, constraining the mass of the se- 
lected candidates to the mass of the B meson. The three- 
body decay kinematics is described by two independent 
di-kaon invariant mass variables (w^m^) = (s23,si3), 
where we order the same-sign kaons such that S23 < S13. 
The signal region contains 1769 B + and 1730 B~ candi- 
dates whose Dalitz plot distribution is shown in Fig. [5] 



III. ISOBAR MODEL FIT 

We perform an extended binned maximum likelihood 
fit to the event distribution in Fig. [5] by binning the 
folded Dalitz plot into 292 non-uniform rectangular bins 
and minimizing the log of the Poisson likelihood ratio, 



Xllr/ 2 = Ei=l Mi 



,SR 



,SR 



In (rt^ R //ii), where 



,SR 



is the number of observed signal region events in the i-th 
bin, assumed to be sampled from a Poisson distribution 
with mean In the limit of large statistics, the Xllr 
function has a \ 2 distribution and can be used to evalu- 
ate the goodness-of-fit. 

The expected number of events in the z-th bin is mod- 
eled as 

m = 2 [e\M\ 2 dm 2 23 dml 3 +R q7! (nf -n^ s )+n^ s , (1) 

J i 

where the first term is the expected signal contribution, 
given by twice the bin integral of the square of the matrix 



FIG. 2: The Dalitz plot of the 1769 B + and 1730 B~ candi- 
dates selected in the signal toes region. The axes are defined 
in the text. 



element multiplied by the signal selection efficiency e, 
determined as a function of (77123,77113) using simulated 
signal events. The integral is multiplied by two because 
we use a folded Dalitz plot. We use the isobar model 
formalism 0, and describe the matrix element M. as 
a sum of coherent contributions, M. = Efc=i-^fc- The 
individual contributions are symmetrized with respect to 
the interchange of same-sign kaons, 1 <-> 2, and are given 



by 



M k 



Pke" 



V2 



(A(s23)P/ fc (cos0 13 ) + {!<-> 2}), (2) 



where p^e % ^ k is a complex-valued decay coefficient, Ak is 
the amplitude describing a K + K~ system in a state with 
angular momentum J and invariant mass y/S23, Pj is the 
Legendre polynomial of order J, and the helicity angle 
6*i3 between the direction of the bachelor recoil kaon 1 
and kaon 3 is measured in the rest frame of kaons 2 and 
3. 

The model includes contributions from the 0(1020) and 
XcO intermediate resonances, which are clearly visible in 
Fig. [3 Following Ref. 01 1 we introduce a broad scalar 
resonance, whose interference with a slowly varying non- 
resonant component is used to describe the rapid decrease 
in event density around m{K + K~) — 1.6GeV/c 2 . Evi- 
dence of a possible resonant S-wave contribution in this 
region has been reported previously |l6l Il7| . however 
its attribution is uncertain: the /o(1370) and /o(1500) 
resonances are known to couple more strongly to 7T7t 
than to KK [l8l|: possible interpretations in terms of 
those states [l!| must account for the fact that no strong 
B ± -> if ± / (1370) or B ± -> AT ± / (1500) signal is ob- 
served in B ± -> K ± 7t ± tt 11 [Hill. The contribution of 
the /o(1710) resonance is included in the fit as a separate 
component and is found to be small. In the following, we 
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designate the broad scalar resonance Xo(1550) and de- 
termine its mass and width directly from the fit. 

The contribution from a spin J resonance with mass 
mo and total width To is described by a relativistic Brcit- 
Wigner amplitude: 



200 



Aj{s) 



rag — s — im (T + Ar(s)) 



(3) 



Fj is the Blatt-Weisskopf centrifugal barrier factor (21 
for angular momentum J: Fq(x) = 1 and Fi(x) e 



x/Vl + x 2 , q h = y/s/A- 



and R represents the effec- 



tive radius of the interaction volume for the resonance; 
we use R = 4.0 GeV" 1 (0.8 fm) In the formulation 

of Eq. (0 , only the centrifugal barrier factor for the de- 
cay of a spin J resonance into two pseudoscalar kaons 
is included; we have ignored the corresponding centrifu- 
gal barrier factor for the two-body decay of a B meson 
into a pseudoscalar kaon and a spin J resonance. The 
effect of this approximation on the parameterization of 
B ± — > K ± 4>(1020), the only component with J > 0, is 
negligible. Unless otherwise specified, all resonance pa- 
rameters are taken from Ref. [lg. The term Ar(s), pa- 
rameterizing the mass dependence of the total width, is 
in general given by Ar(s) = J2i ATj(s), where the sum is 
over all decay modes of the resonance, and Ar^m 2 ,) = 0. 
The Xco has many decay modes, the decay modes of the 
/o(1710) are not well established, and decay modes other 
than K + K~ of the possible Xo(1550) resonance are un- 
known; in all these cases we set Ar(s) = and neglect 
the mass dependence of the total width. For the </>(1020) 
resonance we use Ar(s) = AFi(s) + Ar 2 (s), where 
T t = T B(<j) -> K+K-), T 2 = r o £(0 -» K°K°), and the 
mass dependence of the partial width for the two-body 
vector to pseudoscalar-pseudoscalar decay (/>(1020) — » hh 
is parameterized as 



AIY; 



Ti,: 



qh m$ Fl{q h R) 
qoh \fs Ff(q oh R) 



(4) 
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FIG. 3: The projected m(K + I{-) invariant-mass distribu- 
tions for the best fit: (a) rri23 projection (the inset shows the 
fit projection near the c/>(1020) resonance), (b) mi3 projec- 
tion. The histograms show the result of the fit with BB and 
qq background contributions shown in dark and light gray, 
respectively. 



y / rn|/4~- m" h . 

A large — » i^ ± / (980) signal measured in B ± — > 
J ftT ± 7r ± 7r T [l2l l2fj| . and a recent measurement of gpc/gn, 
the ratio of the /o(980) coupling constants to KK and 
7T7T j23|, motivate us to include an / (980) contribution 
using a coupled-channel amplitude parameterization: 



where q oh 



A 



•/o(980)( 



(5) 



- s - im (g^g n + QkQk) 
where g v = 2/3^/1 - 4mJ±/a + 1/3^1 - 4mys, 



4m,i„/s, and we 



qk — V 2 y i_ ™'ic±/ i ' T1 /'Y*"' 1 '"!r» 
use g K /gir = 4.21 ± 0.25 ± 0.21, m = 0.965 ± 0.008 ± 
0.006 GeV/c 2 and g„ = 0. 165 ± 0.010 ± 0.015 GeV/c 2 ^ . 

We have investigated two theoretical models of the 
nonresonant component |24l l2o| and found that neither 



describes the data adequately. We therefore include an S- 
wave nonresonant component expanded beyond the usual 
constant term as 



M 



NR 



PNR.e 



i<pNR 



V2 



A fit to the ?Tt23 > 2 GeV/c 2 region of the folded Dalitz 
plot of Fig |21 which is dominated by the nonresonant 
component, gives a = 0.140 ± 0.019 GeV"V, = 
— 0.02±0.06 GeV~ 2 c 4 , consistent with no phase variation. 
In the following we fix (3 = 0, and incorporate the A^nr 
contribution over the entire Dalitz plot, thus effectively 
employing the same parameterization as in Ref. M 

We fit for the magnitudes and phases of the decay co- 
efficients, the mass and width of the Xo(1550), and the 
nonresonant component shape parameter a. As the over- 
all complex phase of the isobar model amplitude is arbi- 
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trary, we fix the phase of the nonresonant contribution 
to zero, leaving 14 free parameters in the fit. The num- 
ber of degrees of freedom is 292 — 14 = 278. We perform 
multiple minimizations with different starting points and 
find multiple solutions clustered in pairs, where the so- 
lutions within each pair are very similar, except for the 
magnitude and phase of the Xco decay coefficient. The 
twofold ambiguity arises from the interference between 
the narrow Xco and the nonresonant component, which 
is approximately constant across the resonance. The 
highest-likelihood pair has Xllr = (346.4,352.0); the 
second best pair has Xllr = (362.4,368.7). The least 
significant components are the /o(980) and the / (1710). 
Their omission from the fit model degrades the best fit 
from Xllr = 346.4 to 363.9 and 360.7, respectively. 

The invariant-mass projections of the best fit are 
shown in Fig.|3| The goodness-of-fit is x 2 = 56 for 56 bins 
in the 7/123 projection and x 2 — 66 for 63 bins in the 77113 
projection. The sharp peak in the BB background dis- 
tribution in the TO23 projection is due to the contribution 
from the B ± — > DK^ backgrounds. The fit gives a = 
0.152 ± 0.011 GeV~V, m (X Q ) = 1.539 ± 0.020 GeV/c 2 , 
and T {X ) = 0.257 ± 0.033 GeV/c 2 . The fitted values 
of the shape parameter a and the resonance mass are 
consistent with the values in Ref. ^3] 1 but our preferred 
value for the width is significantly larger. The results 
of the best isobar model fit are summarized in Table [IJ 
where we have also included the results for the XcO com- 
ponent from the second solution in the highest-likelihood 
pair, labeled xlo, and component fit fractions, 



F k = 



J \Mk\ 2 ds2 3 ds 



13 



J\M\ 2 ds 23 ds 

13 



(7) 



where the integrals are taken over the entire Dalitz plot. 
The sum of the component fit fractions is significantly 
larger than one due to large negative interference in the 
scalar sector l2a. 



IV. BRANCHING FRACTIONS AND 
ASYMMETRIES 

To search for possible direct CP violation we extend 
the isobar model formalism by defining charge-dependent 
decay coefficients: 



1±A 



ft e ±iS4> k /2 



(8) 



where Aj~ is the CP asymmetry of the A;-th component, 
and 5(f>k — 4>^ — 0^ ■ We modify the likelihood to be 
the product of the likelihoods for the two charges and 
repeat the fit. The phase of the nonresonant component 
is fixed to zero for both charges. The results are given 
in the last three columns of Table [I] in terms of the fit- 
ted CP asymmetry values, the symmetric 90% confidence 
level intervals around them, and the fitted phase differ- 
ences between the charge-dependent decay coefficients. 



The asymmetry intervals are estimated by fitting Monte 
Carlo simulated samples generated according to the pa- 
rameterized model of the nominal asymmetry fit. There 
is no evidence of statistically significant CP violation for 
any of the components. 

Taking into account the signal Dalitz plot distribution, 
as described by the isobar model fit, the average signal ef- 
ficiency is s — 0.282±0.011, where the uncertainty is eval- 
uated using control data samples, and is primarily due 
to the uncertainties in tracking and particle identifica- 
tion efficiencies. The total branching fraction is £>(-B — > 
K^K^K^) = (35.2 ± 0.9 ± 1.6) x 10~ 6 , where the first 
error is statistical and the second is systematic. The fit 
fraction of the isobar model terms that do not involve the 
XcO resonance is (95.0±0.6± 1.1)% for the best fit, giving 
B(B ± -> i^i^iTF) = (33.5±0.9±1.6)x 10~ 6 if intrin- 
sic charm contributions are excluded. The total asym- 

, • a B(B-^K-K-K+)-B(B+^K + K+K~) 

metry is a C p - B (b-^k- k- k+)+B(b+^k+k+k-) ~ 
(-1.7 ±2.6 ±1.5)%. 

The systematic error for the overall branching frac- 
tion is obtained by combining in quadrature the 3.9% 
efficiency uncertainty, a 1.1% uncertainty on the total 
number of B + B~ pairs, a 0.7% uncertainty due to the 
modeling of BB backgrounds, and a 1.4% uncertainty 
due to the uncertainty arising from the uncertainty on 
the R q q sideband extrapolation coefficient. The 1.5% sys- 
tematic uncertainty for the asymmetry is due to possible 
charge asymmetry in kaon tracking and particle iden- 
tification efficiencies, evaluated using data control sam- 
ples. Where appropriate, the systematic uncertainties 
discussed above have been propagated to estimate the 
uncertainties on the leading isobar model fit results. We 
have also evaluated the systematic uncertainties due to 
the parameterization of resonance lineshapes by varying 
the parameters of all resonances within their respective 
uncertainties. Uncertainties arising from the distortion of 
narrow resonance lineshapes due to finite detector resolu- 
tion and, for candidates containing a 0(1020) resonance 
produced in the qq continuum, due to the kinematic fit, 
have also been studied. 

The values of the partial two-body branching fractions 
are summarized in the fifth column of Tabic [I] Us- 
ing the 6(0(1020) -> K+K-)&nd B(x c0 -> K+K~) 
branching fractions from Ref. we compute S(S ± — > 
#±0(1020)) = (8.4±0.7±0.7±0.1) x 10" 6 and B(B ± -> 
K± Xc0 ) = (1.84 ± 0.32 ± 0.14 ± 0.28) x 10~ 4 , where the 
last error is due to the uncertainty on the 0(1020) and 
Xco branching fractions. Both results are in agreement 
with previous measurements [ToL HH IT2I I2H |2S[ . 

The partial branching fractions for £?± — * #±/o(980) 
measured in the K ± K ± K T and i4' ± 7r ± 7r =F final states are 
related by the ratio 



R = 



g(/o(980) -> K+K-) 

B(/ (980) -» 7T+7T-) 



4 U g-n 



where 3/4 is an isospin factor, and Ik /In is the- ratio 
of the integrals of the square of the /o(980) amplitude 
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TABLE I: The magnitudes and phases of the decay coefficients, fit fractions, two-body branching fractions, CP asymmetries, 
symmetric 90% confidence level CP asymmetry intervals around the nominal value, and the phase differences between the 
charge-dependent decay coefficients for the individual components of the isobar model fit. 

Comp. p <t> (rad) F (%) F x B(B ± —> K^K^K^) A (A min , A ma x)ao% <ty (rad) 

0(1020) 1.66 ± 0.06 2.99 ± 0.20 ± 0.06 11.8 ± 0.9 ± 0.8 (4.14 ± 0.32 ± 0.33) x 10"" 0.00 ± 0.08 ± 0.02 (-0.14,0.14) -0.67 ± 0.28 ± 0.05 

/o(980) 5.2 ± 1.0 0.48 ± 0.16 ± 0.08 19 ± 7 ±4 (6.5 ± 2.5 ± 1.6) x ID -6 -0.31 ± 0.25 ± 0.08 (-0.72,0.12) -0.20 ± 0.16 ± 0.04 

X (1550) 8.2 ± 1.1 1.29 ± 0.10 ± 0.04 121 ± 19 ± 6 (4.3 ± 0.6 ± 0.3) x 10""' -0.04 ± 0.07 ± 0.02 (-0.17,0.09) 0.02 ± 0.15 ± 0.05 

/o(1710) 1.22 ± 0.34 -0.59 ± 0.25 ± 0.11 4.8 ± 2.7 ± 0.8 (1.7 ± 1.0 ± 0.3) x 10 -6 0.0 ± 0.5 ± 0.1 (-0.66,0.74) -0.07 ± 0.38 ± 0.08 

xio 0.437 ± 0.039 -1.02 ± 0.23 ± 0.10 3.1 ± 0.6 ± 0.2 (1.10 ± 0.20 ± 0.09) xKT 6 0.19 ± 0.18 ± 0.05 (-0.09,0.47) 0.7 ± 0.5 ± 0.2 

X ", 0.604 ± 0.034 0.29 ±0.20 6.0 ± 0.7 (2.10 ± 0.24) x 10" 8 -0.03 ± 0.28 - -0.4 ± 1.3 

NR 13.2 ± 1.4 141 ± 16 ± 9 (5.0 ± 0.6 ± 0.4) xlO" 5 0.02 ± 0.08 ± 0.04 (-0.14,0.18) 



given by Eq. JSJ over the B — » KKK and B — > Kirn 
Dalitz plots, and gxlgix is the ratio of the /o(980) cou- 
pling constants to KK and tttt. Using our results and 
those in Ref. [2(j, we get R — 0.69 ± 0.32, where we 
have combined the statistical and systematic errors of 
the two measurements in quadrature. This is consistent 
with R = 0.92 ± 0.07, which we get by evaluating the 
right-hand side of Eq. (jSJ using the values of the ./o (980) 
parameters reported by the BES collaboration |23j . 



V. PARTIAL- WAVE ANALYSES 

We further study the nature of the dominant S-wave 
component by considering the interference between the 
low-mass and the high-mass scattering amplitudes in the 
region m 23 S (1.1, 1.8) GeV/c 2 , m 13 > 2GeV/c 2 . The 
matrix element is modeled as 



M = Pg(f23) e # g ( S23 ) + PNR e - 

V2 \/2 



(10) 



where ps and 4>S refer to the S'-wave and are taken to be 
constant within each bin of the S23 variable and the non- 
resonant amplitude parameterization is taken from the 
fit to the high-mass region. The partial-wave expansion 
truncated at the S-wave describes the data adequately; 
the magnitude of the S'-wave in each bin is readily de- 
termined. Because of the mass dependence of the non- 
resonant component, the phase of the S-wave can also 
be determined, albeit with a sign ambiguity and rather 
large errors for bins with a small number of entries or 
small net variation of the nonresonant component. 

The results are shown in Fig. QJ with the S-wave com- 
ponent of the isobar model fit overlaid for comparison. 
Continuity requirements allow us to identify two pos- 
sible solutions for the phase; the solution labeled by 
black squares is consistent with a rapid counterclock- 
wise motion in the Argand plot around m(K + K~) = 
1.55 GeV/c 2 , which is accommodated in the isobar model 
as the contribution of the X o (1550). 

Isospin symmetry relates the measurements in B^ — > 
i^i^ifT and B° K + K~K° 29]. Our results for the 
K + K~ S-wave can therefore be used to estimate a po- 
tentially significant source of uncertainty in the measure- 
ments of sin2/3 in P° -> 0(lO2O)lf° Hfl HU due to the 
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FIG. 4: The results of the partial-wave analysis of the K + K~ 
S-wave: (a) magnitude squared, (b) phase. The discrete am- 
biguities in the determination of the phase give rise to two 
possible solutions labeled by black and white squares. The 
curves correspond to the S-wave component from the isobar 
model fit. The inset shows the evidence of a threshold en- 
hancement from the fits of the S-wave in the vicinity of the 
K + K~ threshold and in the region around the <^(1020) reso- 
nance. 



contribution of a CP-even S-wave amplitude. We per- 
form a partial-wave analysis in the region ni23(K + K~) G 
(1.013, 1.027) GeV/c 2 , which we assume to be dominated 
by the low-mass P-wave, due to the contribution of the 
</>(1020) resonance, and a low-mass S-wave. The matrix 
element is modeled as 



M= ps + Mpl elMs23)cos 



(11) 



where the low-mass S-wave is taken to be constant over 
the small S23 interval considered. The fit results for the 
P-wave are shown in Fig.[3 with a Breit-Wigner fit of the 
0(1020) resonance overlaid for comparison. For the S- 
wave we get p| = (3.4± 2.5) x 10 2 GeV _4 c 8 and compute 
its fraction in this region using Eq. J7J) to be (9 ± 6)%. 

We also consider the region 2m K + < m(K + K~) < 
1.006 GeV/c 2 , in the immediate vicinity of the K + K~ 
threshold. The contribution of the 0(1020) resonance 
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FIG. 5: The results of the partial-wave analysis in the 0(1020) 
region for the P-wave: (a) magnitude squared, (b) phase. The 
discrete ambiguities in the determination of the phase give rise 
to two possible solutions labeled by black and white squares. 
The curve corresponds to a Breit-Wigner fit of the 0(1020) 
resonance. 



tail in this region is suppressed by the centrifugal barrier 
and is estimated to be smaller than 10%. We fit p| = 
(6.1 ± 1.6) x 10 2 GeV~ 4 c 8 for the magnitude of the S- 
wave in this region. The fits in the vicinity of the K + K~ 
threshold and in the region around the 0(1020) resonance 
indicate a threshold enhancement of the S'-wavc, which 
is accommodated in the isobar model by the contribution 
of the /o(980) resonance as shown in the inset of Fig. ^ 



VI. CONCLUSIONS 

In conclusion, we have measured the total branching 
fraction and the CP asymmetry in £?± — > K^K^K^ . An 



isobar model Dalitz plot fit and a partial-wave analysis 
of the K + K~ S-wave show evidence of large contribu- 
tions from a broad A"o(1550) scalar resonance, a mass- 
dependent nonresonant component, and an /o(980) reso- 
nance. The ratio of B ± -> X ± / o (980) two-body branch- 
ing fractions measured by BABAR in £>± — > X ± 7r ± 7r =F 
and -B± — > K^K^RT is consistent with the measure- 
ment of QkIQt: by the BES collaboration, albeit with 
large errors. Our isobar model fit results are substan- 
tially different from those obtained in Ref. 01 due to 
the larger fitted width of the Xq(1550) and the inclusion 
of the /o(980) component in the isobar model. Our re- 
sults for the B(B ± -> #±0(1020)) and B{B ± -> #±x c o) 
branching fractions are in agreement with the previous 
results from BABAR [ToL which they supersede, and 
from other experiments |l2l |28| . We have measured 
the CP asymmetries and the phase differences between 
the charge-dependent decay coefficients for the individual 
components of the isobar model and found no evidence 
of direct CP violation. 

We wish to dedicate this paper to Prof. Richard E. 
(Dick) Dalitz, inventor of the "Phase Space Plot" (as he 
called it) with which his name is so intimately linked, in 
tribute as much to his qualities as a gentleman as to his 
gifts as a physicist. 
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1 



TABLE I: Fit fractions matrix of the best fit. The diagonal elements Fkk correspond to component fit fractions shown 
in the paper in Table I. The off-diagonal elements give the fit fractions of the interference terms defined as Fu = 
25R / MkMtds 2 3ds 13 / J\M\ 2 ds 2[i ds 13 . 



Fm x 100% <p /q(980) Xq(1550) /o(1710) Xco NR 

<f> 11.8 ± 0.9 ± 0.8 -0.94 ± 0.18 ± 0.11 -1.71 ± 0.36 ± 0.24 0.01 ± 0.10 ± 0.03 0.11 ± 0.02 ± 0.05 3.54 ± 0.38 ± 0.40 

/o(980) 19 ± 7 ± 4 53 ± 12 ± 7 -4.5 ± 2.9 ± 1.2 -0.9 ± 0.2 ± 0.5 -85 ± 21 ± 14 

X (1550) 121 ± 19 ± 6 -30 ± 11 ± 4 -1.1 ± 0.3 ± 0.5 -140 ± 26 ± 7 

/ (1710) 4.8 ± 2.7 ± 0.8 -0.10 ± 0.07 ± 0.07 4 ± 6 ± 3 

Xco 3.1 ± 0.6 ± 0.2 3.9 ± 0.4 ± 1.9 

NR 141 ± 16 ± 9 



TABLE II: Fit parameters correlation matrix of the best fit. 







/o(980) 


Ao(1550) 


/o(1710) 


XcO 


NR 




p <i> 


P <t> 


p <j> m T 


p 4> 


P <f> 


p a 




1 2 


3 4 


5 6 7 8 


9 10 


11 12 


13 14 


1 

2 


1.00 -0.10 
-0.10 1.00 


-0.14 0.01 
0.69 0.61 


0.00 -0.00 0.02 -0.00 
0.38 -0.42 -0.64 0.46 


-0.00 0.01 
0.03 -0.20 


-0.01 0.01 
0.02 -0.29 


-0.09 -0.07 
0.34 0.21 


3 
4 


-0.14 0.69 
0.01 0.61 


1.00 0.72 
0.72 1.00 


0.33 -0.63 -0.44 0.37 
0.68 -0.80 -0.50 0.67 


0.05 -0.37 
0.19 -0.51 


0.05 -0.30 
0.07 -0.41 


0.73 0.49 
0.70 0.56 


5 
6 
7 
8 


0.00 0.38 
-0.00 -0.42 

0.02 -0.64 
-0.00 0.46 


0.33 0.68 
-0.63 -0.80 
-0.44 -0.50 

0.37 0.67 


1.00 -0.52 -0.14 0.94 
-0.52 1.00 0.39 -0.63 
-0.14 0.39 1.00 -0.30 

0.94 -0.63 -0.30 1.00 


0.45 -0.59 
0.28 0.75 
0.25 0.05 
0.17 -0.63 


0.05 -0.42 
-0.02 0.30 
-0.01 0.14 

0.04 -0.39 


0.46 0.41 
-0.52 -0.29 
-0.07 0.00 

0.34 0.26 


9 
10 


-0.00 0.03 
0.01 -0.20 


0.05 0.19 
-0.37 -0.51 


0.45 0.28 0.25 0.17 
-0.59 0.75 0.05 -0.63 


1.00 0.04 
0.04 1.00 


0.05 -0.19 
0.00 0.22 


0.36 0.42 
-0.39 -0.20 


11 
12 


-0.01 0.02 
0.01 -0.29 


0.05 0.07 
-0.30 -0.41 


0.05 -0.02 -0.01 0.04 
-0.42 0.30 0.14 -0.39 


0.05 0.00 
-0.19 0.22 


1.00 0.40 
0.40 1.00 


0.10 0.17 
-0.35 -0.30 


13 
14 


-0.09 0.34 
-0.07 0.21 


0.73 0.70 
0.49 0.56 


0.46 -0.52 -0.07 0.34 
0.41 -0.29 0.00 0.26 


0.36 -0.39 
0.42 -0.20 


0.10 -0.35 
0.17 -0.30 


1.00 0.92 
0.92 1.00 



TABLE III: Phase differences matrix of the best fit. The matrix elements are denned as 4>ki = 4>i — <t>k- The last row corresponds 
to the results of the best fit shown in the paper in Table I. The other phase differences are derived from the fit results. 



4>kl 



/o(980) 



Xp(1550) 



/o(1710) 



XcO 



NR 



4> 

/o(980) 
X (1550) 
/o(1710) 

XcO 

NR 





2.51 ± 0.17 ± 0.11 
1.70 ± 0.27 ± 0.07 
-2.71 ± 0.35 ± 0.13 
-2.27 ± 0.35 ± 0.12 
2.99 ± 0.20 ± 0.06 



-2.51 ± 0.17 ± 0.11 


-0.81 ± 0.25 ± 0.12 
1.07 ± 0.36 ± 0.17 
1.50 ± 0.33 ± 0.17 
0.48 ± 0.16 ± 0.08 



-1.70 ± 0.27 ± 0.07 2.71 ± 0.35 ± 0.13 2.27 ± 0.35 ± 0.12 

0.81 ± 0.25 ± 0.12 -1.07 ± 0.36 ± 0.17 -1.50 ± 0.33 ± 0.17 

-1.88 ± 0.19 ± 0.09 -2.32 ± 0.22 ± 0.09 

1.88 ± 0.19 ± 0.09 -0.44 ± 0.30 ± 0.13 

2.32 ± 0.22 ± 0.09 0.44 ± 0.30 ± 0.13 

1.29 ± 0.10 ± 0.04 -0.59 ± 0.25 ± 0.11 -1.02 ± 0.23 ± 0.10 



-2.99 ± 0.20 ± 0.06 
-0.48 ± 0.16 ± 0.08 
-1.29 ± 0.10 ± 0.04 
0.59 ± 0.25 ± 0.11 
1.02 ± 0.23 ± 0.10 




